Averaging rheological quantities in descriptions of soft glassy materials 
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Many mean-field models have been introduced to describe the mechanical behavior of glassy 
materials. They often rely on averages performed over distributions of elements or states. We here 
underline that averaging is a more intricate procedure in mechanics than in more classical situations 
such as phase transitions in magnetic systems. This leads us to modify the predictions of the recently 
proposed SGR model for soft glassy materials, for which we suggest that the viscosity should diverge 
at the glass transition temperature T g with an exponential form r\ ~ exp( T J 4 T )■ 

PACS: 62.20.Fe, 64.70.Pf, 83.80.Hj 
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I. INTRODUCTION 

Concentrated colloidal suspensions of soft particles 
(emulsions, slurries) tend to display glassy behavior when 
concentrated beyond a threshold Jl|-^1. Despite their 
great physical and chemical diversity, these Soft Glassy 
Materials often share many mechanical features. This 
has prompted theoreticians to construct generic models 
[p|-pO| , independent of the microscopic features of the sys- 
tems, and often borrowed from statistical physics. The 
outcome of such a procedure is often a description in 
terms of a large distribution of relaxation processes act- 
ing simultaneously, a picture supported by experimental 
facts for many glassy systems ]ll],[l2| . 

Any macroscopic mechanical response of the system 
then results from an averaging over these distributions. 
In this paper we stress the importance of the averag- 
ing procedure when one deals with mechanics. Aver- 
aging procedures have already been discussed for non- 
mechanical situations in the non-ergodic glassy phase 
(e.g. statistical mechanics of spin glasses [0). How- 
ever averaging of mechanical fields is a priori a different 
story (l3) . This also affects the description of the ergodic 
fluid phase (e.g. the divergence of the viscosity when one 
approaches the "glass transition" from the fluid side) . 

To make our point we focus here on the recently in- 
troduced SGR model (Soft Glassy Rheology) || which 
incorporates mechanics in a simple picture of the glass 
transition introduced by Bouchaud. The system is taken 
as a collection of blocks, that evolve according to stochas- 
tic equations, driven by the imposed shear rate. A control 
parameter of the model is the effective temperature that 
we will denote here T. A transition occurs at T g below 
which the system is no more ergodic (and the fluid dis- 
plays a yield stress). A remarkable feature of the model 
is that the viscosity diverges at 2T g (not T g ), so that in 
the interval T g < T < 2T g the material is a power-law 
fluid. This large domain of viscosity divergence, in the 
ergodic phase, is a feature absent in most other models 



for soft glassy materials [^pJio[] . 

After general comments about averaging procedures in 
mechanics (section II) , we show below that this feature is 
a consequence of the adopted averaging procedure (sec- 
tion III). We propose an improvement of this procedure, 
that leads to a Vogel-Fulcher divergence of the viscosity 
at T g : r\ ~ exp( T ^ T ). We end by a short discussion and 
by drawing a few research directions (section IV). 



II. VARIOUS AVERAGES FOR 
RHEOLOGY/MECHANICS 

Let us start by clarifying the crucial role of the aver- 
aging procedure in rheology or mechanics, very different 
from the one used in many classical areas, such as mag- 
netic systems, where the imposed field H is well-defined 
, and the average magnetization is simply the algebraic 
mean of the local magnetizations. In rheological prob- 
lems, the fields are imposed through the boundaries of 
the system, so that often the shear rate and the shear 
stress play more or less symmetric roles, which results in 
a wide variety of situations. 

To clarify the message, let us focus on the simple case 
of a system built of elements (blocks) of high viscosity 
?7i and elements (blocks) of low viscosity 772 (Figs, la 
and lb), and estimate the average viscosity of the system 
through different averaging procedures. 

A first average corresponds to the situation where the 
same shear rate 7 is imposed to all elements, and one 
averages the resulting stresses, which leads to an average 
viscosity 



iv)* = + 0- - <Pi)v2 



(i) 



where <f>\ is the proportion of blocks of high viscos- 
ity. Strictly speaking this procedure describes blocks ar- 
ranged in layers disposed in parallel (Fig. lc). A second 
kind of average corresponds to submitting all blocks to 
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a sparse dispersion of highly viscous blocks in a less 
viscous matrix (Figl.b), it yields the Einstein formula 
{i]) OP ~ 772(1 + §</>i), and thus similar physics to the 
constant stress picture (77) a ~ 772(1 + (j>i) (Fig. Id). 
For a sparse dispersion of fluid inclusions in a highly 
viscous matrix (Fig. la), (ri) OP ~ 771 (1 — §(^2), a re- 
sult qualitatively akin to that of the constant shear rate 
picture (77}^ = 771(1 — (fe) (Fig. lc). Of course the 
Olroyd-Palierne formula is unable to describe accurately 
the properties of intermediate mixtures (roughly <fi in the 
range 30 — 70%) as it neglects spatial structures that 
could result in percolation effects of great consequence 
on the mechanical behavior. 



FIG. 1. Schematic pictures of : (a) a dispersion of fluid 
inclusions 2 in a highly viscous matrix 1, (b) a dispersion of 
highly viscous inclusions 1 in a fluid matrix 2, (c) a geome- 
try corresponding to averaging at fixed strain rate 7, which 
leads to a highly viscous behavior as in (a), (d) a geometry 
corresponding to averaging at fixed stress with a qualitative 
behavior similar to (b). 



III. CONSEQUENCES FOR THE SGR MODEL 

Implicitly the SGR model of Sollich et al. [|| considers 
an ensemble of elements of viscosity distributed along: 



r)(E) = 77 exp(£;/r) 



(5) 



the same stress <j, and averaging the shear rates. This 
results in a different average viscosity 



(v) a = OMr 1 + (i - Mife 1 ) -1 



(2) 



which describes layers in series (Fig. Id). 

The behavior according to the first kind of average is 
dominated by the highly viscous blocks: the average vis- 
cosity (77)^ diverges if 771 diverges, whatever the relative 
volume fractions. This is obviously a reasonable repre- 
sentation of a highly viscous medium with low-viscosity 
inclusions (</>i large, Fig. la), but a poor one for a system 
of highly viscous inclusions dispersed in a fluid matrix of 
low viscosity 772 (Fig. lb). On the contrary, the latter 
case (<pi small, Fig. lb) is better represented by a fixed 
stress average (Fig. Id), as (77) is dominated by the low 
viscosity blocks. 

Clearly none of these two simple averages properly 
describes all situations. Consequently effective medium 
theories have tried to provide better alternatives. Among 
those, let us quote the Olroyd-Palierne formula [[uj for a 
matrix of viscosity rj m with various kind of dilute spher- 
ical inclusions i (volume fraction fa, viscosity 77,). This 
formula states that the effective viscosity of the resulting 
composite fluid is (neglecting surface tension effects): 
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with 



OP 



H, = 



1 - T,ifaHi 



2rj t - 2rj„ 
2r\i + 377,, 



(3) 



(4) 



The above formula provides a convenient interpolation 
between the two limiting physics mentioned above. For 



where E characterizes the activation barrier limiting the 
yield of these blocks. The distribution of these barriers 
is modeled by n(E) ~ exp(—E/T g ), which naturally in- 
troduces T g as a (Glass) transition temperature. Indeed, 
for T > T g , the states are populated according to Boltz- 
mann statistics at steady state p(E) ~ n(E) exp(E/T) 

E) , which ceases to 



which gives p(E) = T T 9 exp( 
be normalizable as T — > T+. 

The average viscosity used in this model is the constant 
rate one, which leads for T > T g to : 



T-T g 



dEp(E)r](E) = / dE 



,T — Tn 



S.E 



TT„ 



(6) 



which diverges for T ->• 2T+ (i.e. not for T -> T+). With 
this averaging procedure, the viscosity is thus infinite in 
the range T g <T < 2T g , above the glass transition tem- 
perature. The corresponding mechanical behavior was 
found to be that of a power-law fluid with a temperature 



dependent exponent {a) ~ 7 T n . 

Let us now explain why, on physical grounds, we ex- 
pect this prediction to be incorrect, and the viscosity 
of the system to remain finite for T g < T < 2T g . For 
example, focus on the 10% most viscous blocks (i.e. 
those of viscosity larger than r] c = 770 exp(E c /T) with 



E r = - 



TgT 

T—T„ 



log(O.l)). These few very viscous blocks 
are embedded in a matrix of viscosity smaller than rj c , as 
all its blocks have a viscosity weaker than r\ c . The effec- 
tive viscosity of the whole fluid, composed of this matrix 
with sparse quasi-solid inclusions, is thus smaller than 
T/c times a geometric factor related to the 10% volume 
fraction (e.g. 77 < ^f§- q r\c in the Olroyd-Palierne model), 
and the actual position of the inclusions. The viscosity 
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is thus finite for T g < T < 2T g and can not diverge faster 
than T] c = % exp(-^r- log(10)) as T -> T+ . 

Switching to the other limit and using a fixed stress 
average, leads to averaging the inverse of the viscosities 
(or the inverse of the relaxation times). The continuous 
version of Equation (2) yields an average viscosity: 



T 



T —T„ 



(7) 



which displays a weak (algebraic) divergence at T g 

We now propose a physically sounder recipe based on 
the Olroyd-Palierne formula and a self-consistent approx- 
imation. We take a representative but small proportion 
4> of the blocks as "inclusions" , while the rest is the "ma- 
trix". As the matrix and the whole system are similar, 
formula (3) gives an average viscosity {rj) OPSS — r\ m : 
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2n(E) + 3 (rj) 



(8) 



(9) 
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The average viscosity (rj} p SS is thus implicitly given by 
Y,i<f>iHi = (a similar approximation was used in another 
context JTJ): 
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(10) 
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Although the actual numerical factors (2 and 3) may be 
discussed, we believe that this formula gives a physically 
correct average, intermediate between the two limiting 
cases mentioned above: (77)^ > (r])oPSS — ( 7 ?)<r' an< ^ 
with the appropriate behavior in the two limits. 

Applying this to the SGR model of Sollich et al. leads 
to a divergence of the viscosity at T g . Indeed if we 
write the average viscosity at temperature T in the form 



i^OPSS 
for E* is: 



t?7o exp(i?*/T), then the implicit equation 



dE p(E) 



exp^^ 



3/5 



(11) 



Given that J °° dEp(E) = 1, this suggests J^l dEp(E) ~ 
2/5 and thus exp(-.E*^f) ~ 2/5. This corre- 
sponds to a viscosity diverging for T — > T+ as ~ 
log(5/2)). This divergence, obtained by 



3770 exp( 5 r^7 

the above rough analysis, is confirmed by numerical in 
spection of (10) which gives 



with a prefactor ~ 0.72 instead of 2/3. This exponen- 
tial divergence is consistent with our qualitative analysis 
below equation (6). 

Of course the 5/2 is not to be taken too seriously, but 
we think that a Vogel-Fulcher like divergence of the vis- 
cosity at T g is a sound result, within the hypothesis of 
the SGR model. 



IV. DISCUSSION 

In conclusion, great care should be taken when choos- 
ing an averaging procedure for the mechanical properties 
of a system close to a "glass transition" . Fixed stress and 
fixed shear rate averages are only extremes that limit a 
large set of possibilities. A physically appealing inter- 
mediate procedure for elements of distributed properties, 
adapted from effective medium theory, has been proposed 
in equation (10). 

Focusing on the SGR model, we have shown that us- 
ing different averages seriously modifies the law charac- 
terizing the divergence of the viscosity. Using the recipe 
of equation (10) we have proposed that a Vogel-Fulcher 
divergence could quite naturally arise. Revisiting the 
model, to incorporate the effects of the inhomogeneity 
of the driving strain rate field 7 on the dynamics, is a 
direction that we wish to explore. 

It would also be interesting to analyze thoroughly what 
kind of averages are implicitly made within the dynamical 
rules used in other models for such systems 0,^,|lO|. For 
example the simple models of |?],|| predict an algebraic 
divergence of the viscosity. A natural question is whether 
this is the result of an implicit averaging of the inverse 
of the relaxation times as in equation (7). 

We have considered here the simplest (in principle) sit- 
uation of an ergodic system (i.e. above the effective glass 
transition temperature), where neglecting spatial struc- 
tures and adopting mean-field approaches is most likely 
to apply. Analysis of non-ergodic situations is left for fur- 
ther work. Eventually let us clarify that we have focused 
here on the sole viscosity to make our point clear, but 
that a more complete description of rheological proper- 
ties would have to deal with effective theories for assem- 
blies of visco-elastic (and possibly plastic) elements [[l4| . 
Such effective theories most likely will have to deal with 
the importance of spatial correlations and with the role 
of percolating structures. 

Acknowledgments: we thank S. Obukhov for fruitful 
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